Dynamic radiation force of acoustic waves on solid elastic spheres 



Glauber T. SilvH 

Centro de Pesquisa em Matemdtica Computational, 
Universidade Federal de Alagoas, Maceio, AL, Brasil, 57072-970 
(Dated: February 2, 2008) 

The present study concerns the dynamic radiation force on solid elastic spheres exerted by a 
plane wave with two frequencies (bichromatic wave) considering the nonlinearity of the fluid. Our 
approach is based on solving the wave scattering for the sphere in the quasilinear approximation 
within the preshock wave range. The dynamic radiation force is then obtained by integrating 
the component of the momentum flux tensor at the difference of the primary frequencies over the 
boundary of the sphere. Results reveal that effects of the nonlinearity of the fluid plays a major role 
in dynamic radiation force leading it to a parametric amplification regime. The developed theory 
is used to calculate the dynamic radiation force on three different solid spheres (aluminium, silver, 
and tungsten). Resonances are observed in the spectrum of the force on the spheres. They have 
larger amplitude and better shape than resonances present in static radiation force. 

PACS numbers: 43.25.+y 



I. INTRODUCTION 

Among nonlinear effects of acoustic wave propagation 
in fluids is worth noticeable radiation pressure (force) 
and streaming 0. These phenomena come together as 
an incident wave encounters an obstacle. They are phys- 
ically different but may combine producing an effective 
force upon the obstacle. Usually these phenomena can- 
not be separated in actual fluids ;3j. When thermovis- 
cous effects are negligible, streaming can be ruled out. 
With respect to time dependence, radiation force is ei- 
ther static or dynamic. Static radiation force is a time 
averaged quantity produced on a target by a monochro- 
matic acoustic wave. The dynamic counter-part is an 
oscillatory force exerted on a target by a bichromatic 
acoustic wave. Here we consider that this force oscillates 
in time at the difference frequency of the wave, i.e., the 
difference between the primary wave frequencies. 

The utilization of dynamic radiation force in science 
and technology has broadened. It has been applied for 
measuring ultrasound power of ultrasonic transducers M , 
inducing oscillation in bubbles and liquid drops 0, 
and exciting modes in capillary bridges |7|. Dynamic 
radiation force has become the underlying principle in 
some elastography techniques such as shear wave elas- 
ticity imaging 8] and vibro-acoustography 9J. Using vi- 
brometry, the viscoelasticity of transparent materials can 
be determined by locally measuring the vibration of a em- 
bedded sphere induced by dynamic radiation force [Icj . 

In a pioneer work, King ^1 calculated the static radi- 
ation force exerted on a rigid sphere by a plane traveling 
wave in a lossless fluid. Study on the static force of com- 
pressible and elastic spheres followed up [l2L IT^ | . It was 
promptly noticed that the static radiation force exhibits 
fluctuations related to the mechanical properties of the 
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sphere. Effects of viscosity and thermal conductivity of 
the host fluid on the static radiation force over a rigid 
sphere were also theoretically accounted 0] . For elastic 
solid media, a theoretical formulation of radiation force 
problem can be found in Ref . |15| . 

The influence of the nonlinearity of the fluid on static 
radiation force is negligible at least for lossless fluids 
where the nonlinear parameter B/A < 100 [T3.IT7|. The 
dynamic radiation force exerted on a solid sphere by a 
plane traveling wave has been theoretically investigated 
in Ref. [l^. Experimental validation of the theory is 
found in Ref. 19]. In both works, the difference fre- 
quency was assumed to be very narrow. On that account 
the magnitude of dynamic and static radiation forces are 
alike. Dynamic radiation force exerted on an object is 
directly caused by the pressure wave at the difference fre- 
quency acting on the object. It is known that the wave 
at the difference frequency might be subjected to para- 
metric amplification [20j . This suggests that dynamic ra- 
diation force itself should achieve a regime of parametric 
amplification. Silva et al. |2l| verified that the dynamic 
radiation force on an acrylic sphere can be greatly am- 
plified as the difference frequency is increased. 

The aim of this work is to analyze the dynamic radi- 
ation force exerted on a solid elastic sphere in the para- 
metric regime. This analysis should be carried consid- 
ering the nonlinearity of the fluid in the problem. To 
the best of the author's knowledge this problem has not 
yet been treated. To accomplish this task we have to 
solve the acoustic scattering for the sphere in at least the 
quasilinear approximation. In what follows we present, 
based on regular perturbation technique, the wave equa- 
tions in the quasilinear approximation for lossless fluids. 
The solution for a bichromatic plane traveling wave is 
derived from these equations within the preshock wave 
range. Using the Mie scattering theory we obtain a de- 
scription for the scattered waves by the elastic sphere. 
The yielded result is used in an vector integral formula 
obtained for the dynamic radiation force. We apply the 
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developed theory to analyze the dynamic radiation force 
on three different elastic spheres, namely, aluminium, sil- 
ver, and tungsten. A connection of dynamic radiation 
force and the Resonance Scattering Theory (RST) is es- 
tablished. The spectrum of the dynamic radiation force 
for the spheres is numerically evaluated. Resonances are 
observed in the spectrum of the force. They depend in 
a unique manner on the mechanical and elastic parame- 
ters of the sphere. Thus, the resonances in the spectrum 
could be used for characterization of materials. Further- 
more, their magnitude are up to 81 dB higher than those 
found in static radiation force. 



nates (r, 9, tp) is given by [21 



II. NONLINEAR WAVE PROPAGATION AND 
SCATTERING 



Consider a homogeneous and isotropic fluid with adi- 
abatic speed of sound Co, in which thermal conductivity 
and viscosity are neglected. The fluid has infinite extent 
and is characterized by the following acoustic fields: pres- 
sure p, density p, and particle velocity v = — \/<p. The 
function <f> is the velocity potential and V is the gradient 
operator. These fields are function of the position vector 
r and time t. At rest, these quantities assume constant 
values p — po, p = po, and v = 0. The acoustic fields are 
governed by the dynamic equations of ideal fluids. By us- 
ing the regular perturbation technique, one can expand 
velocity potential in terms of the Mach number e <C 1 
as <j) — £<f>^ + e 2 (f>^ + 0(e 3 ), where the super-indexes 
indicate the order of the potential functions. Any analy- 
sis of radiation force has to be done considering at least 
the quasilinear approximation, i.e., second-order acoustic 
fields. 

To obtain the dynamic radiation force over an object, 
we have to solve the scattering of a bichromatic wave for 
the object. Consider a bichromatic traveling plane wave 
of finite-amplitude with primary frequencies lo\ and u>2 
propagating along the z-axis. The wave is formed by a 
sinusoidal excitation at z = 0. Due to nonlinear nature 
of wave propagation, a third wave at the difference fre- 
quency UJ21 = — w\ arises in the fluid. Suppose the 
bichromatic wave hits a sphere of radius a placed in the 
z-axis at the distance zo from the acoustic source. The 
elastic sphere has density pi, comprcssional and shear 
speed of sound denoted, respectively, by c c and c s . To 
simplify our analysis we consider k2izg 3> 1 and a <C zq. 
Let us denote the complex amplitudes of the potential 
functions (incident plus scattered) at the frequencies uji , 

LU2, and u>2i by, respectively, <fii, (f>2, and <p2i- The total 
potential for scattering problem can be found using the 
Mie scattering theory. Accordingly, the velocity poten- 
tials in terms of spherical waves in the spherical coordi- 
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where A m = c G e~ lk ^ / \2k m ) , A 21 = 1Co z e- lk ^ z « /4, i 
is the imaginary unit with k m = Lu m /co and &21 = LO21/CQ. 
The quantity 7 = l+B/A, where B/A being the so-called 
nonlinear parameter of the fluid. The functions hn^ 1 and 
hh, are, respectively, the first- and second-kind spherical 
Hankel functions of nth-order, while P n is the Legendre 
polynomial of nth-order. 

The quantity S n is the modal scattering function. To 
find this function we apply the solid elastic-fluid bound- 
ary conditions on the sphere surface. These conditions 
yield a linear system whose solution is 
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The matrix elements of this equation are found in 
Ref. [2I. 



III. THEORY OF DYNAMIC RADIATION 
FORCE 

Let the Fourier transform of a function g(t) be denoted 
by F\g\- The inverse transform is represented by J-^ 1 . 
The dynamic radiation force exerted by a bichromatic 
wave on an object is given by [l8| 
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The quantity pov^ v^ 1 ) is the Reynolds' stress tensor, 
while p^ is the second-order pressure given by 



P {2) = e 2 po 



1 /<90« 
2? \ dt 
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We restrict our analysis to waves in which U21 < u>i, 
otherwise we may have superposition of the dynamic ra- 
diation force at W21 and other forces caused by the waves 
at the primary frequencies. 
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The second integral in the right-hand side of Eq. (JSJ 
depends on the variation of the object boundary caused 
by the incident wave. The evaluation of this integral for 
a moving object is a hard task because this is a self- 
coupled problem. We have to calculate the force on the 
object in first-order by integrating the linear pressure on 
the boundary of the object. However, the pressure may 
deform or move the boundary of the object. We shall 
show an approximate calculation of this term for a rigid 
spherical target. 

Consider a rigid sphere of radius a and density pi sus- 
pended in a fluid. The origin of the coordinate system 
coincide to the center of the sphere at rest. We assume 
the sphere oscillates only along the z-axis. The variation 
Az in the first-order approximation of the normal vector 
of the sphere is given by 

n = e, Aze z + 0{e 2 e z ) 1 

where e r is the radial unit vector and e 2 is the unit vec- 
tor along the z direction. We have Az = Szie~ tuJlt + 
Sz 2 e~ ZUJ2t , Szi and Sz2 are the amplitudes of oscillation 
of the sphere. The second integral of Eq. (JSJ) might be 
approximated to 



where p„ = epodc^m' /dt is the linear pressure and 
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The amplitudes of oscillation are given in terms of the 
mechanical impedances of the sphere Z\ and Z 2 at the 
frequencies u>t and ui 2 as [23| 



■ 4 3 

1-7TO W„ 

* > 



Po- 



IX" 



pi 



1,2. 



The quantity x m — w m a/co is the size factor of the 
sphere. Equations JJJ and (jSJ) give us an approxima- 
tion to calculate the first term of the amplitudes of the 
dynamic radiation force in Eq. 

The dynamic radiation force is calculated as follows. 
With the obtained first- and second-order potential fields 
in Eqs. Q and J2J), the first term in the right-hand side 
of Eq. (JSJ is calculated performing the surface integral. 
The second integral of this equation can be evaluated 
using Eqs. Q), and 10. Accordingly, the dynamic 
radiation force on the sphere is given by 
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where E — s 2 poCq/2 is the energy density at the acous- 
tic source and Y21 is the radiation force function. This 
function might be decomposed as Y 2 i — Y a + Y b + Y c , 
where 
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where R n {x) = i n [h (2) (x) + S^h^ (x)} and R' n (x) = 

i n [hn^ (x) + S n {x)h^ (x)]. The symbol ' means the 
derivative of the function with respect to its argument. 
The function Y21 depends on the size factors x 2 and x\ 
with the restriction x 2 \ < x\ < x 2 . 

The function Y a is proportional to the inverse of the 
radiation impedance of the sphere at the fundamental fre- 
quencies oj\ and lo 2 . This term is divergent when the size 
factors xi and x 2 approach to zero. In the case x 1; x 2 > 1 
we may disregard the contribution of Y a to the radiation 



force function. Furthermore, lim Y a = 0. The contribu- 

£21^0 

tion of function Y\, brings up the regime of parametric am- 
plification in the dynamic radiation force. It has a domi- 
nant role whenever the size factor x 2 \ is not a small quan- 
tity. It also depends linearly on the nonlinear parame- 
ter 7 and the ratio zq/cl. This function has the follow- 
ing asymptotic behavior \Yb\ ~ 2zoja~ 1 (l + x^}\ . 

The amplification of the radiation force with zq does 
not increase indefinitely. On that account two aspects 
should be asserted. The distance z should be within the 
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preshock wave range otherwise Eqs. (JTJ and J5J are no 
longer valid. Diffraction and fluid viscosity may alter the 
dependence of radiation force with the distance Zq (sec 
Ref. The function Y c has a prominent contribu- 

tion to the dynamic radiation force whenever x mn <C 1. 
When X21 — the incident wave becomes monochromatic 
and the dynamic radiation force in Eq. 1)1011 reduces to 
its static counter-part as obtained in Ref. 

Hereafter, we focus on the situation of which X21 is 
not a very small quantity and x\, X2 > 1. Therefore, the 
radiation force function becomes Y21 = Yb- In that case, 
the integral formula of dynamic radiation force in Eq. (J5J 
is reduced to 

f 2 i(*) = - [[ P 2 2 i(t)ndS, (14) 

J J S 

where p^i = £ 2 Pod4> 2 2 i /dt with cf>$ being the second- 
order velocity potential at the difference frequency. No- 
tice the case for which X21 <C 1 has been analyzed else- 
where [3. 



IV. RESONANCE SCATTERING THEORY AND 
DYNAMIC RADIATION FORCE 

In our previous analysis, we have found that the dy- 
namic radiation force function reduces to Y21 = Yj,. Ac- 
cording to the Resonance Scattering Theory (RST), the 
response of an elastic sphere to a plane wave excitation 
can be decomposed into a modal rigid body background 
interfering with resonance contributions |22j. The reso- 
nances appear due to the eigenvibration modes of some 
types of surface waves (Rayleigh and whispering gallery 
waves) present on the sphere as the result of the plane 
wave excitation. 

To isolate the resonance contributions to the dynamic 
radiation force, we follow the approach developed in 
Ref. (2j|. Accordingly, the resonance scattering function 
is given by 



where S„\x) = ' (x) / ' (x) is the modal scattering 
function of a rigid sphere. We thus define the resonance 
radiation force function of an elastic sphere as 

F 2 ( r S) =X2lS[ IeS \x 2 l)h^(x2l), (16) 

for which we dropped the factor 72:0/0 because it is not 
relevant to resonance analysis. The resonance radiation 
force function is proportional to the backscattered wave 
(6 = 7r) of mode n = 1 (the dipole mode). Further- 
more, this function carries an unique information of the 
sphere in terms of its mechanical parameters, including 
the elastic constants of the material. 



V. NUMERICAL RESULTS 

We evaluate the radiation force function for three dif- 
ferent solid spheres, which are made aluminium, silver, 
and tungsten. The physical parameters under consid- 
eration of these materials are summarized in Table [I] 
The fluid medium is water whose parameters are pa — 
1000 Kg/m 3 , c = 1500 m/s, and 7 = 6. 

TABLE I: Physical parameters for the spheres after Ref. |26|. 



Speed of sound 



Material Density Compressional Shear 





[Kg/m 3 ] 


[m/s] 


[m/s] 


Aluminium (Al) 


2700 


6374 


3111 


Silver (Ag) 


10500 


3704 


1698 


Tungsten (W) 


19250 


5221 


2887 



In Fig.p we plot the spectrum of the function Y21 — 
a (7- z o) _1 ^2i for the three spheres. This function varies 
within the range < X21 < 30 incremented by steps of 
0.01. We also consider x\,X2 > 1, therefore, the contri- 
butions from Y a and Y c to the radiation force function are 
negligible. The region < X21 < 5 corresponds to the 
crossover between the static radiation force (X21 = 0) 
and the fully developed dynamic radiation force 221 > 5. 
To see how larger is the magnitude of the dynamic radi- 
ation force in this region compared to that of the static 
counter-part consider z^ja = 1000. The magnitude of 
the static radiation force function for all spheres is about 
the unit j2(J, while that for the dynamic radiation force 
function in X21 > 5 is around 12000. Therefore, the mag- 
nitude of the fully developed dynamic radiation force is 
81 dB higher than that of the static radiation force. The 
spectrum of dynamic radiation force function follows the 
behavior of the radiation force function for a rigid sphere 
in which c c , c s , p\ — > +00, except for the the appearance 
of well defined resonances. The labels Al-n, Ag-n, and 
W-n indicate the nth resonance in the radiation force 
function for, respectively, the aluminium, the silver, and 
the tungsten spheres. Only the first three resonances are 
highlighted in the graphs. 

The spectrum of the resonance radiation force func- 
tion for the aluminium sphere is plotted in Fig. [21a. We 
plot the phase of the corresponding resonance scattering 
function in Fig. [21b. Whenever the spectrum undergoes 
through a resonance or an antiresonance (destructive in- 
terference of two adjacent peaks (2?|]) the phase shifts by 
7r. The almost constant behavior of the phase near to 
the first peak reveals that the peak is not related to a 
resonance of the sphere. 

In Fig. [SJa, we have the spectrum of the resonance 
radiation force function for the silver sphere. The phase 
of the resonance scattering function is shown in Fig. 0b. 
We have a similar behavior of the spectrum as observed 
for the aluminium sphere. 
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FIG. 1: Magnitude of the radiation force function Yi\. The 
labels Al-n, Ag-n, and W-n indicate the nth resonance in 
the radiation force function, (a) Aluminium, (b) Silver, (c) 
Tungsten. 
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FIG. 3: The spectrum of the resonance radiation force func- 
tion of the silver sphere, (a) Magnitude, (b) Phase of the 
resonance scattering function S^i^ ■ 
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FIG. 2: The spectrum of the resonance radiation force of the 
aluminium sphere, (a) Magnitude, (b) Phase of the resonance 
scattering function S^i^ ■ 



FIG. 4: The spectrum of the resonance radiation force func- 
tion of the tungsten sphere, (a) Magnitude, (b) Phase of the 
resonance scattering function S^i^ ■ 



The spectrum of the resonance radiation force func- 
tion for the tungsten sphere is presented in Fig. 0] The 

behavior of the resonances as well as the phase of s[ lc ^ 
follow those of the aluminium and the silver spheres. 



VI. DISCUSSION AND CONCLUSION 

We have carried out a theoretical analysis of dynamic 
radiation force caused by bichromatic waves on a target 
taking into account the nonlincarty of the fluid medium. 
In the analysis, the object target must be within the 
preshock wave range. To obtain the force exerted on an 
solid elastic sphere we had to solve the scattering problem 



in the quasilinear approximation using the partial wave 
expansions given in Eqs. Q and The incident and 
scattered fields were used in Eq. JSJ) yielding the radia- 
tion force. It was demonstrated that dynamic radiation 
force is caused by parametric amplification of the differ- 
ence frequency wave. This is experimentally supported 
by Ref . [21| . Further, the force is proportional to the non- 
linear parameter of the fluid 7. We therefore presented 
a much simpler formula for dynamic radiation force in 
Eq. IH|). 

The obtained result of the dynamic radiation function 
in Eq. I|12|l permitted a straight connection with this 
force and the Resonance Scattering Theory (RST). To 
the best of the author's knowledge such connection has 
not yet been explored before, neither theoretically nor 
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experimentally. The spectrum of the dynamic radiation 
force as a function X21 exhibited well defined resonances, 
which have larger amplitude (about 81 dB higher) and 
are better shaped than those present in static radiation 
force |26(. The resonances are caused by eigenvibrations 
in some surface waves on the sphere at the difference 
frequency W2i- The RST allowed us to separate the res- 
onance part of the radiation force function as given in 
Eq. (|16|l . We, thus, analyze the resonance spectrum for 
three different spheres made of aluminium, silver, and 
tungsten. Furthermore, the study of dynamic radiation 
force on spherical shells can be readily done by using the 
present theory. The location, amplitude, and width of 
these resonances depend upon the mechanical and the 
elastic properties of the sphere. Hence, the spectrum 
information could be used as the material signature in 
a noncontact acoustic spectroscopy technique. In this 
method, the induced vibration by the dynamic radiation 
force on a specimen could be measured by a laser vi- 
brometer yielding the spectrum of the force [l^. From 
the spectrum it is possible to extract the elastic proper- 
ties of the specimen in a similar manner as done in the 



Resonant Ultrasound Spectroscopy psj . 

We have not considered thermoviscous effects in the 
fluid. By neglecting such effects we ruled out acoustic 
streaming in the boundary of the sphere. Usually stream- 
ing has the tendency of altering the total effective force 
acting on a target. Moreover, the dependency of dynamic 
radiation force with temperature was not accomplished 
here. Variations in the temperature may change the lo- 
cation and width of the resonances in the spectrum of 
the dynamic radiation force. 

In conclusion, this paper presented a comprehensive 
theoretical analysis of the dynamic radiation force on 
solid spheres. The spectrum of this force carries use- 
ful information of the sphere related to its elastic and 
mechanical properties. 
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